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(54) Method and apparatus for calculating color gamuts 



(57) A method and an apparatus are disclosed to 
obtain a gamut description of a multidimensional color 
reproduction device. The method is based on the trans- 
formation of several types of colorant boundaries that 



map to the boundaries of the color gamut of the device 
in color space. 
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Field of the invention 

The present invention relates to devices and methods for image rendering by means of multidimensional color 
reproduction devices. 

Background of the invention. 

.^^^^^^^^^•^^ 

coSs for a™ primer. The variables for the colorants are denoted as o,. * c n wth n the d.mens.on of the 

is complete* characterized by its coiorant gamut with a number of colorant limitations and the printer 

See of this close relationship between an n-ink process and the printer m«xlel^e ^ 
model are also defined for the n-ink process. With inverting an n-.nk process « meant that the corresponding pr.nter 

S^sformSn of an n-ink process to co.or space on the other hand is equhralent to the transformation of the cor- 

n-h* oraS S setting n-m colorants to constant values in the printer model of the n-.nk process, the cctorarrt hmrta- 
Sons a elrSmX this m-ink process if the colorant gamut of the nvink process is restricted by colorant I m.te^ ons 
<o obSn2 by seSg the n-m colorants in the colorant limitations of the colorant gamut of the n-,nk process to the,r con- 
stant values. orocess with m < n means that n-m out of the n colorants in the printer 

Zn^ueTthey can have according to the colorant gamut of the n-ink process and the , «, lorant Mons are .nher- 
45 ited bv the m-ink process. The m-ink process is called the extracted m-.nk process of the n-mk process. 

K mink processes extracted from an n-ink process, and the n-m out of the n colorants in the pnnter jnodel of the n- 

nTpScess^e Replaced by their minimum or maximum value as defined by the colorarj ^gamut of the n-,nk process. 

he m-ink process is called an m-ink boundary process of the n-ink process. In general there are 
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A ctorS is completely determined by its boundaries. In a 3-dimensionaJ color space these boundaries are sur- 
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1 faces, and hence there will be surfaces in colorant space that map to the color gamut boundaries. In general several 
kinds of colorant surfaces should be taken into account but for well-behaved printers the color gamut is, by definition, 
exactly determined by the transformation of the colorant domains of the 2-ink boundary processes of the n-ink process 
to color space. It can be shown that this condition is equivalent with saying that in a n-ink process, a colorant combina- 

b tion lies inside the color gamut if there are at least three colorants that can be changed both positively as negatively. 
Such surfaces in colorant space are called physical colorant boundaries that are transformed by the printer model to 
the physical color boundaries in color space. The envelope of the physical color boundaries determines the boundaries 
of the physical color gamut of the n-ink process. In the case of an s-dimenslonal color space, the dimension of the 
boundaries is s-1 . Physical colorant boundaries correspond in this case to the domains of the (s-1 ) -ink boundary proc- 

10 esses of the n-ink process. The physical color boundaries are obtained by transforming the domains of the (s-1) -ink 
boundary processes of the n-ink process to color space. 

The color gamut of a printer is certainly one of its most important characteristics. An explicit representation of this gamut 
is highly desirable, for example, to enable the optimal reproduction of images that contain colors that fall outside the 
color gamut of the printer. 

is Because our invention is based on the use of a printer model, a few more words of explanation are given on the subject 
of printer models. 

The Neugebauer equations 

20 Different mathematical expressions have been proposed to mode! printing processes, along which the Neugebauer 
equations are certainly the most known ones. 
) The Neugebauer equations reflect how color originates in halftone printing. The printing with three inks c 1 , eg and C3 by 
means of three halftone screens results theoretically in eight possible combinations of colorant overlap, called the Neu- 
gebauer primaries. 

25 If the size of the halftone dots is small enough not to be resolvable by the eye, the resulting color is a weighted sum of 
the tristimulus values of these primaries. For a 3-ink process, the Neugebauer equation for the X tristimulus value is 
given by 

^(C, , C 2 , C 3 ) =a VV^lV* 3 1 ^1 +a 2^2 +a 3^3 +a 12^12 +a 13^13^ a 23^23 +a 123^123 

30 

If it is assumed that the relative positions of the halftone dots are random, the Neugebauer equations can be calculated 
from the Demichel equations that predict the fraction of each combination of the three inks as a function of their respec- 
tive dot percentages c 1( c^ and C3, and this leads to the Neugebauer equations in their most often encountered form: 

a i=(c 1 )(1-c 2 )(1-c 3 ) 
a 2=( 1 - c i)(c 2 )0-C3) 

40 

) a 3=0-C 1 )(1-C 2 )(C 3 ) 

a 23=(1-C 1 )(C 2 )(C 3 ) 

45 a 13=(C 1 )0-C 2 )(C 3 ) 

a 12=(Cl)(C 2 )(1-C 3 ) 
a 123=(C 1 )(C 2 )(C 3 ) 

50 

in which c-j , c 2 and C3 represent the dot percentages of the three inks. The equations for the Y and Z values are obtained 
by replacing the X values by the corresponding Y and Z values respectively. The extension of the Neugebauer equa- 
tions for n inks is straightforward. 

According to one interpretation, the Neugebauer equations predict the X, Y and Z tristimulus values by means of trilin- 
55 ear interpolation from the corresponding tristimulus values of the 8 Neugebauer "primaries". If the Neugebauer equa- 
tions are used as multilinear interpolation formulae, they can be used to model any color device in any color space. 
Working out the Demichel coefficients in the Neugebauer equations and rearranging the terms results in the following 
set of polynomial expressions : 
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X = k 0+ k 1 C 1+ k 2 C 2+ k 3 C 3+ k 12 C 1 C 2+ k 13 C 1 C 3+ k 23 C 2 C3 + k 123 C 1 C 2 C3 

Y = 1 0+ l ,c 1+ l 2 c 2+ l 3 c 3+ l 12 o n c 2+ l n3 c ,03+1230203+1 n23 o 1 o 2 c 3 
Z = m 0 + mi o 1 +m202+m 3 03+m 1 20 1 02+m 13 o 1 o 3 +m 23 o 2 o 3 +m 123 o 1 o 2 o 3 

^sleadstoaninterp^ 

be calculated from the X. Y and Z values of the ^^K gir eq^ can be graphically rep- 
££*n technique on a large set a curved surface in the col- 

Improving the precision of the Neugebauer printing model. 

ThTwemraalo" « •» Newt*"" ■*"•»»» "^T^^nSheSicsplof -pttowlMtrllinMrlnw- 

25 ability of Color, pp. 32-37. 

Calculation of the color gamut boundaries. 
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Caicuianon ot xn» wiwi 

Anexam,eofafirstheuris«cme^^ 

crevices", presented at TAGA 94 by »^^ lC ^S^ col0 rant oombinations are evaluated by 
According to this method, the colorant gamut J, for example, different C1E lightness val- 

means of aprinter model. The resutong .^^"^eTe^an approximation of the shape of the color gamut 
ues. The envelope of the colors that ,,e w *' n on * 

at the corresponding constant ■^^^"JS s SSd •> achieve a reasonable description of the color 
the first place a large set of samples m f "" 8 ^ J^sing time for this method increases exponentially 
gamut boundaries. Hence this method is ^ slow and ' ^J?, a m^men^tically point of view it is difficult to determine 
with the number of colorants. A second P roW ^« th c t\SCTseTof discrete points. Particularly the detection of 
exactly the boundaries of the color gan^ 

holes or several contours per lightness level is rafter conpUcMeo. n t te too large so that 
rcbustbecause it may faiito detect the to S^^c\ *S are not seen as connected by the 

neighboring points in colorant space are transformed to colors .n co. h» 

the color gamut and making the color increment smaller, * he ^^ ^ ^ oM and a sufficient number of 
be determined with any desired precision. By selecbnga , BuH«.enUj ™* esfimated This method, however, 

directions for the color change, lindanes ^^spaoe so the printer model has to be inverted 
also has a number of drawbacks. ^^"SJ^^Sonrt combinations give rise to the same ^color 
Till now this is almost always done by iterative methods. conclusions will be drawn for the 

only one solution will be found. If M ^^^? e ^^Sdomafn. A second problem is that it can 
color gamut K at least one of the other '^^^^^ cot « am have to be determined. A third draw- 
fail when the color gamut is concave or w^P* 0 ^ ^^culations since many more points need to be 
back is that the heuristic nature of the procedure leads to very a 



4 



EP0 763 930 A1 



1 evaluated than are obtained for the description of the gamut boundary Itself. A faster version of this method is repre- 
sented in the article "The color gamut obtainable by the combination of subtractive color dyes", Journal of Imaging sci- 
ence, Vol. 30, No. 1, 1986 by N. Otha. Per lightness plane a border point is searched for and the next border point is 
detected by looking radially around this point From this new border point the next border point is detected again by radi- 
s ally looking around it. This last procedure is repeated until the first border point is reached again. In this way concave 
contours can be detected, but nevertheless this iterative method is still quite slow, the printer model has to be inverted 
and multiple contours are not detected. 

A third approach is described in the article by Masao Into, "Fast Algorithm for Computing Color Gamuts", Color 
Research and Application, Vol. 18, No. 5, Oct. 1993. In order to calculate a gamut descriptor of a 3-ink process, thefol- 
10 lowing procedure is performed: 

Step 1 : The lightness values at the eight vertices of the dye delimited solid are calculated for zero and maximum 
amounts; 

Step 2: For a given lightness, the locations of several equilightness points on the edges of the dye delimited solid 
is are computed; 

Step 3: Loci connecting these points on the faces of the dye delimited solid are computed; 
Step 4: The loci are mapped into color space to produce the boundary of the color gamut. 

According to the article, the loci of the equilightness points on the edges of the dye solid are computed as follows. Since 
20 at the edge of the color solid one of the colorants is equal to its minimum or maximum value, there are effectively only 
two colorants available as variables to achieve the desired lightness value. One of these variables can be used as a 
) parameter, leaving only one value left to iterate on. One dimensional Newton-Raphson is used to determine the value 
of this remaining variable that results in the desired lightness value. 

The above method is claimed to work very fast but its weak point is that it relies on Newton Raphson iteration. Depend- 
25 ing on the printer model, there may be several possible solutions, but only one solution will be found with iterative meth- 
ods. Nevertheless, if several start points would be used to find ail the solutions, it is not possible to know which colorant 
combination belongs to (which one of) the previous solution(s), i.e. no information is available about the connectivity 
between the different solutions. 

Moreover, for unconventional processes some boundaries of the color gamut are not obtained by mapping the edges 
so of the colorant cube to color space. An example of such a 3-ink process was found in printing yellow (y), cyan (c) and 
green (g) inks. The XY2 and corresponding CIELAB values of the Neugebauer primaries of this process are given in 
Table 1 . The first column corresponds to the 8 possible combinations of ink overlap. An evaluation reveals that for this 
process approximately half of the realizable colorant combinations yields colors that lie outside the color gamut as it is 
defined by the mapping of the edges of the colorant cube to color space. Figure 2 shows a cross section of the color 
35 gamut in CIELAB at a constant L* value. The solid lines correspond to the intersection with the mapping of the edges 
of the colorant cube to color space. The shaded area represents colors that lie outside the gamut defined by the phys- 
ical boundaries. The outer curve represents the intersection with the object color solid, also known as the MacAdam 
limits. 
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L* 


• 

a 




white 


100 • 0 


1 00 0 


100.0 


100.0 


0.0 


0.0 


y 




66 2 


10.5 


85.1 


-3.1 


80.1 


g 


10 7 


29 .8 


7.0 


61.5 


-43.1 


51.2 


c 


25.7 


33.6 


69.2 


64.7 


-29.7 


-37.8 


yg 


14.1 


19.6 


1.7 


51.4 


-30.0 


64.6 


yc 


11.1 


19.8 


6.9 


51.6 


-51.0 


34.4 


8 C 


4.7 


11.2 


4.6 


39.9 


-60.9 






| 3.0 


J 6.8 


1.1 


31.4 


-49.3 





Table i : 

non- conventional 3 -ink process. 



To WM the «*» 0»™t « "J"*** «* "J^15^ eited problems rsmsin. Ho«evsr. 

with moTe than three colorants this is certainly no longer the case. 
Objects of the invention. 

ical colorant boundaries. 

» •. so owed of M rnv**>n to <Wo*. «U th. 00.0™* boundaries «», wnstom K ft. Moodoriss of » 00I0, 
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gamut for a n-ink process with n larger than the dimension of the color space, taking into account several colorant 
limitations. 

Further objects and advantages of the invention will become apparent from the description hereinafter. 
Summary of the invention 

The above mentioned objects are realised by the specific features according to claim 1. Preferred embodiments of 
the invention are disclosed in the dependent claims. The method according to the current invention is especially well 
suited for use in a three dimensional color space, also referred to as tristimulus color space. If colors are described in 
a four dimensional color space or a color space with an even higher dimension, the method according to the current 
invention is preferentially applied by reducing the number of dimensions to three This may be accomplished by not 
using one or more dimensions, or by using the sum of two or more variables in color space, or any other function applied 
to such variables, such that finally three useful independent variables remain. 

In a three dimensional color space, a color gamut is a solid. For describing this solid, the color gamut is preferentially 
intersected by a plurality of surfaces in color space. A surface is a set of points according to two degrees of freedom, 
and may be curved. A plane is a specific embodiment of a surface. A more accurate description of the color gamut may 
be obtained by intersecting it by more surfaces. For each surface, all points on the surface, which belong to the color 
gamut, are situated in one or more connected regions. Each such region is enclosed by a closed contour. A contour is 
a curve in color space, having one degree of freedom. In the most simple case, all points of the color gamut, situated 
on one surface, may be enclosed by just one closed contour. By specifying contours, each contour being situated in one 
surface, enclosing colors belonging to the color gamut of the reproduction device, the entire color gamut may be 
described. By knowing these contours, it may be easily decided whether a specific color belongs to the color gamut or 
not 

Knowledge of this is important in color reproduction, in order to decide how colors outside the color gamut will be repro- 
duced. 

In a preferred embodiment, each surface is planar. In a specific embodiment, parallel planes in color space, each hav- 
ing a constant lightness may be used to intersect the color gamut In that case, each closed contour is a planar closed 
curve, situated in a plane parallel to each other closed contour. 

According to the current invention, a plurality of closed contours In color space may be constructed. We have found that 
it is more efficient to find curves in colorant space, corresponding to the closed contours in color space. Each curve is 
preferentially based on a sub-ink process, in which three inks or colorants are selected. Each curve is then transformed 
from colorant space to color space. Such a transformed curve gives a portion or segment of a closed contour in color 
space, or even the full closed contour. According to the method of the current invention, one or more contour segments 
may be concatenated. If "one" contour segment is "concatenated"* this means that its starting point is connected with 
its endpoint, forming a closed contour. If two contour segments are concatenated, this means that the endpoint of the 
first contour segment is connected with the starting point of the second contour segment and the endpoint of the second 
contour segment is connected with the starting point of the first contour segment Concatenation is preferentially done 
in colorant space, but it may equivalents be done in color space, after transforming the curve to color space. 
In the broadest sense, Neugebauer cells are solids, as a portion from a colorant space, such that each point belonging 
to the colorant gamut belongs to at least one Neugebauer cell. Such Neugebauer cells may thus have common points, 
regions or solids. 

In a specific embodiment, the Neugebauer cells partition the colorant gamut in mutually disjunctive portions. This 
means that the intersection between each two Neugebauer cells is empty or comprises at most one surface, in a more 
preferred embodiment, each Neugebauer cell is a cuboid, i.e. a rectangular parallelepiped. A specific embodiment of a 
cuboid is a cube. The colorant gamut may be split cubes, all having the same size. 

in the broadest sense, "Neugebauer equations" as referred to in the current invention are equations having at least one 
non-linear term. A non-linear term is e.g. a quadratic term, cubic, etc. or a cross-product between two or more linear 
variables (e.g. Ci*c£. AH sets of equations, all having constant or linear terms only, are excluded from the concept of 
"Neugebauer equations". 

Preferentially, Neugebauer equations are polynomials that express color values in terms of bilinear, triiinear or quadratic 
colorant values. In a specific embodiment, Neugebauer equations are limited to the formulae for X, Y and Z given under 
the title The Neugebauer equations" in the background of the invention, which are referred to as triiinear equations. 
By a 2-colorant space is meant the domain of a 2-ink process, in the invention QV95149, the color gamut was obtained 
by transforming the physical colorant boundaries to color space. For well-behaved processes, such as the cmyk offset 
printing press, the correct color gamut is obtained, but for some unconventional processes, certain colorant combina- 
tions may transform to colors in color space that fall outside the physical color gamut. In this invention all the other kinds 
of colorant boundaries are defined that map to possible boundaries of the color gamut 

Suppose that colors are represented in a 3-dimensional space. For 3-ink processes, only one other kind of boundaries 
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20 Fig. 4 
Fig. 5 
Fig. 6 



> *~i boundaries In the case of more than three colorants, also a hybrid form of 
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sented. 

Detailed description of the invention. 

£nd vertical lines are the a* and b* axes intA . inC0 , 0 r space. The incremental colorant changes 

B migrate in to pia-tenaHted to an, wo * llhjW ^^^^^^weNeutiebaue, primates" 
T«> natural a»mdanwm«»r*4^«^ 

by the appending claims. 

Color gamut description hodsh to 

case the color gamut is described with contours. 
Color gamut boundaries of a 3-ink process 

the color gamut boundaries in color space. 
Physical boundaries 

^ on. ^ «* - « r ssr-is cJSEttsEsssK- 

Soa IW a 3-ir* proeaaa. to gamut In ar»ra™ W»« Zmu^^Sl conatent * altor <ha minimum or maxrmum 
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' It is possible to indicate that, assuming that a color can be rendered by only one set of colorants, the physical bounda- 
ries indeed correspond to boundaries in color space. This is demonstrated by means of Figure 3 and 4. For a point A 
inside the colorant gamut, there exists a corresponding point A 1 in color space. The Neugebauer equations are allowed 
to be linearized for incremental changes of the inks, and hence does a small change dc-i of the ink c-i correspond to a 
5 change in color space that can be characterized by a vector V-, . Similarly does a change of the other two inks and C3 
result in color changes V 2 and V 3 . Change of an ink q in the opposite direction results in a color change opposite to the 
corresponding vector Vj. 

If, however, the point A lies on a boundary of the colorant gamut, as in Figure 4 where A lies in the colorant plane with 
c 1 = c imax • ti* en ,he corresponding color A* can only be changed in the plane defined by the vectors V 2 and V 3 , and 
10 for negative changes of V 1 . 

As a result, only migrations to half of the color space can be achieved by altering colorants and from this follows that 
the point A' is indeed a point of the color gamut boundary. 

If the Neugebauer equations are used as a printing model, it is shown in the invention GV951 49 that the physical bound- 
aries have the shape of quadric surfaces in color space. The cross section of such surfaces with any other plane yields 
is a conic section and this is in the invention QV95149 the basis of the analytical method to calculate gamut descriptors 
consisting of, for example, slices with constant Y values. 



Natural boundaries 



20 For well-behaved three ink processes, such as the cmy offset printing press, the boundaries of the color gamut are 
indeed obtained by transforming the six boundary planes of the colorant gamut to color space, and hence the color 
gamut boundaries do correspond to the physical boundaries. 

It was found, however, that for certain unconventional printing processes, colorant combinations can be found inside the 
colorant gamut that are yet transformed into colors that fall outside the physical color gamut. The latter can only mean 
25 that in such cases the "real" color gamut has to be larger than the physical color gamut, and that additional boundary 
surfaces have to be looted for that define the color gamut and th&t go through the colorant gamut. These additional 
boundary surfaces are introduced here and are called natural boundaries, as they relate more to the way that color is 
created in the printing process than to limitations on the realizable range of the colorants. 

The derivation of the natural boundaries is as follows. If the color A' in Figure 3 lies inside the color gamut, then this 
30 means that an incremental color change can be achieved in any direction by means of an appropriate selection of three 
colorant changes (dci,dC2,dC3). Inversely, if the color A* lies on the boundary of the color gamut, this means, by defini- 
tion, that migrations from this color are not possible in any direction. A first example of such a situation are the physical 
boundaries. A second possible situation is depicted in Figure 4: because the three vectors V 1b V 2 and V 3 happen to lie 
in a plane, only incremental color migrations within this plane are possible, and hence is the point A' part of a natural 
35 boundary. The surface itself is found by working out the condition that the three vectors \A| , V 2 and V 3 are to lie in a 
plane, or - which is equivalent - that they be linearly dependent 

The XYZ components of the color change V A caused by the colorant deviation d<H are expressed by: 

V J dX ^ BY dz\ 
1 V ' dc x ' dc x ) 



45 and similarly are the XYZ components of the color changes V 2 and V 3 given by: 



50 



1/ J J*X. dz \ 

V * [dc* dc 2 ' dcj 
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r \dc 3 ' 3c 3 ' dc 3 ) 



ing determinant be set to zero: 



into requiring that the foHow- 
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dX 


dY 


dc x 




dX 


ay 


dc 2 


dc 2 


dX 


dY 


k 3x 3 


3c 3 



=0 
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quadratic function in C3 : 

(mo+mi c l+mi1 c 2 n+ m 112 c 2 ^m^m^c* 2+ m 122 c n c 2+ m ^ = 0 
Tnisindicatesttat^^ 

form two natural boundary surfaces. ^ tn6 Neugebauer primaries are given injablel. 

Raure 5 shows two natural boundary surfaces of the ^^^J^Zamcesscan be larger than what is predicted by 
nmaylookSradlctory at first that the ^f^^^SS^SL lies inside the coior gamut. This 
Sysical boundaries, and hence, that a color on « P^^^^Sorirt combination that renders that same color 
oaradox can only be explained by the presence of at M«" ^""T^ colorant cnang es can be found indeed 
STes not lie on a physical boundary nor on a natural bound** JJ^K^, ^ , imite set by the first colorant 
Z result in color changes in any <£*£"**2; exists'between the presence of mulfcple col- 

u^of a cotor ol particular rrrtereat Tire teal Ihrea ^"^^'^TTT,^, ^ colorant cornblnatlona thai each 
Ocelot AftS. -TM»r, apac*. ^ " " "^'^^Sries In Flovra 2. C*r T also Ilea 

r^ra^-r". 6 ^^^^ — - 

expected this color can be separated into 1wo sets of not color "5" lies outside the physical boundary yet 

course lies on the physical ^^^^Z^^ofaAora^. Final* the.astcolorl.es on 

mthe above example, a physical boundary -sfound J^^J^ ^ coior gamut. In general, the cdor 
lope of both of them. 
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Table 2 : Dot percentages of the ink process in Table 1 that 



result in the corresponding colors. 



Color gamut boundaries of 4-ink processes 
Natural boundaries 

For a 4-ink process, an incremental change of the ink q results in a color change represented by the vector V ; in color 
space. A first kind of boundary originates if the 4 vectors Vj corresponding to the 4 inks lie in one plane. This leads to 
the extension of the definition of a natural boundary as it was introduced for a 3-ink process. Mathematically, a natural 
boundary is encountered if the rank of the following matrix is lower than 3 

f dX/dc x dx/dc 2 dx/dcs dx/dc^ 
dY/dc x dY/dc 2 dY/dc 3 dY/dc c 
pZ/dc^ dZ/dc 2 dZ/dc 3 dZ/da^ 



This is equivalent to 
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det 



and 



{dX/dc x dX/dc 2 3X/6c 3 
dY/dc^ dY/dc 2 dY/dc 3 
K dZ/dc r dZ/dc 2 dZ/dc 3 ) 



=0 



10 



15 



det 



ax/ac, ax/ac 2 ax/ac^ 

dY/dc x dY/dc 2 dY/dCt 
$Z/dc x 3Z/dc 2 dZ/dCt) 



=0 



25 



30 



^rZrS. ^r? A and the second matrix contains terms with 

K 2 c 2 1 +L 2 c 1 +M 2 = 0 
The solutions of both equations are given by : 



with Kl K 2> L,, I* Mi. M 2 polynomials in c 2 . 



c i= 



35 



40 



E^hothsolu^ 

? y suSace of the 4- ink process. Evaluation of a vanety J four vec tors in a plane is even less likely 

in practical situations. This is not ^ ""^^^."S^ re safe to state that for the practical calcu- 
tooccurLntheco^ 
lation of the color gamut of 4-ink processes me naxu 

Hybrid boundaries 

» ™i«t that lies inside the 4-dimensional colorant 
From the previous paragraph it is clear that for practical, purposes, a point that has 

• , ^ mM m. n8 lonal.Fora4.lnkp«>c^«>«^*s™ su * eo,0,a ' ,, ' 5 ° unaan ' S 
Because Sy color M kfe *n Die color ^ 0< " e °^ ^^toM seated W on th B0 =^bo»d»»s of 

affected by the physical limitation on the fourth ink and nence ao w 
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Physical boundaries 

in the case of a physical boundary of one of the eight 3-ink boundary processes, the color gamut limitation is only due 
to physical restrictions on inks, including on the fourth one. The physical boundaries of a 3-ink process are 2-ink bound- 
s ary processes, so the corresponding physical boundaries of the 4-ink process are also 2-dimensiona! with the third and 
fourth ink set at their minimum or maximum value. 

For "well behaved" processes, such as the cmyk process, the color gamut is not affected by natural boundary condi- 
tions. In that case, the gamut is only defined by the physical boundaries, i.e. 24 different 2-ink boundary processes. 

10 Color gamut boundaries of n-ink processes 

in the case of an n-ink process with n>4, also natural and hybrid boundary types have to be considered. However, in 
practice it is safe to take into account only the physical and natural boundaries of the 3-ink boundary processes. This 
means that the color gamut of an n-ink process is the union of the color gamuts of all its 3-ink boundary processes. 
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Color gamut descriptor 



As the color gamut of an n-ink process is the union of the color gamuts of its 3-ink boundary processes, a color gamut 
descriptor for 3 colorants has to be presented. For a 3-ink process, only two kinds of boundary surfaces have to be 
so taken into account, i.e. the physical and natural colorant boundaries. If a gamut descriptor consisting of a set of contours 
is used, for example the cross sections of the color gamut with equiluminance planes, only a method of the intersection 
with the natural planes has to be presented. The intersections of the physical boundaries with equiluminance planes is 
presented in the invention GV95149. 

26 Neugebauer model 

If the printer is modeled with the Neugebauer equations, the natural boundary planes are described by 



Ac 2 3 + Bc 3 + C = 0 



with 

A = ko+k^+^cfe+k^c^ 



B = l 0 +llC 1 +i 11 C 2 1 +l 1n2 C 2 1 C2+l2C2+ll22ClC 2 2+ll2ClC 2 C3 

35 C = mo+m-j &j +m«j «j c^ +m<j -j 2^1 C2+m 2 C2+m 2 2C^2 + n r h 22^1 c^ 2 +n"H 2C1 Cj> 

Hence there are two solutions for c 3> i.e. As a result there will be two surfaces in colorant space, one for 



„ -S±//g 2 -4AC 
C 3 = ZA 



the +sign and one for the -sign, i.e. the two natural sign boundaries. Per natural sign boundary, the c^ domain is sam- 
pled and the corresponding value for C3 is calculated. If the c 1 C2 domain is sampled per percentage, the c-, C2 combi- 

45 nations can be represented with a gray value image with 101 rows and 101 columns. 

The gray value is equal to the corresponding C3 value, in a second step the pixel values are replaced by the luminance 
value. This value is obtained by evaluating the Y-Neugebauer equation for the colorant values corresponding to the a, 
C2 position and the corresponding c 3 value. The intersection of the natural boundary plane with a constant luminance 
plane is obtained by thresholding the image with the given luminance value and detecting the contours by making use 

so of a contour follower. This results in contours in the 3-dimensional colorant space, that are transformed into the lumi- 
nance plane by making use of the Neugebauer equations. The envelope of these contours and the contours obtained 
by intersecting the physical color boundaries with the luminance plane gives the intersection of the color gamut of the 
3-ink process for the given luminance plane. 

55 Localized Neugebauer model 

The gamut of a 3-ink process modeled with the localized Neugebauer model is obtained by taking the union of the gam- 
uts of all the Neugebauer cells. Because each Neugebauer cell is modeled with the Neugebauer equations for three 
colorants, the gamut calculation in the previous paragraph can be used to determine its color gamut. 
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General printer mode! 

^ Nfiuaebauer equations with any accuracy if a 

equations. 

Co.orgamutdescnptorwHhco.o^Wtations recalcul atedandthecoloraat«imita- 
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Limitation on the number of simultaneously printed colorants 

The color gamut of processes in which a color can only be obtained with a limited number of colorants can be deter- 
mined by making use of the procedures presented in the invention QV95149. 

In most practical situations, the number of simultaneously printed colorants is" limited: A typical example is hifi-color 
printing with 7 inks, it is almost impossible to construct a model for 7 inks because too many colorant combinations have 
to be printed and measured. Therefore the printer model of an n-ink process consists of a number of submodels, that 
each represent the printer for a limited number of colorants. In the case of hifi-color printing, the process could be mod- 
eled with a number of 4-ink processes, and as a result the color gamut is the union of the color gamuts of these 4-ink 
processes. In general, if an n-ink process is modeled with a number of m-ink processes with m < n, the color gamut is 
the union of the color gamut of the m-ink processes. If there are colorant limitations for the n-ink process, they are inher- 
ited by the m-ink processes. 

If however, a model is available for the n-ink process, and only m inks may be printed simultaneously, only m-ink bound- 
ary processes of the n-ink model are taken for which the remaining h-m inks are equal to 0%. The color gamut of the n- 
ink process is the union of the color gamuts of these m-ink boundary processes. 

Having described in detail preferred embodiments of the current invention, it will now be apparent to those skilled 
in the art that numerous modifications can be made therein without departing from the scope of the invention as defined 
in the following claims. 

Claims 

1 . Method for obtaining a description of a color gamut for a color reproduction device, describing said color gamut by 
a plurality of closed contours, each situated on one surface iri color space, comprising the following steps ; 

- characterizing said color reproduction device by a printer model for transforming a set of colorant values in col- 
orant space to a set of color values in color space ; 

extracting at least one 3-ink process from said printer model ; 

- per extracted 3-ink process, determining at least one sufface in colorant space by setting the determinant of 
the Jacobian of the printer model of said extracted 3-irik process to zero ; 

- per extracted 3-ink process determining a plurality of curves in colorant space, each curve being based upon : 

said printer model ; 

- one surface in color space ; and, 

- at least one colorant boundary or at least one colorant limitation or at least said surface in said colorant 
space ; 

transforming each curve to a contour segment in color space, according to said printer model ; 

per extracted 3-ink process taking the envelope of all the contour segments in color space to form at least one 

closed contour ; 

- taking the union of said closed contours from extracted 3-ink processes. 

2. Method according to claim 1 , wherein said curve is determined by evaluation of a function in colorant space or by 
a method comprising the steps of : 

establishing a plurality of samples in a 2-colorant space ; 

- computing from each sample at least one value according to said printer model ; 

- finding said curves in colorant space by following contour lines having the same level for said computed value. 

3. Method according to any of claims 1 to 2, wherein transforming said curve comprises the steps of : 

- sampling said curve to obtain a plurality of samples iri colorant space; 

- transforming each sample to color space, accordihg to said printer model. 

4. Method according to any of claims 1 to 3, wherein the colorartt values are limited to a colorant gamut and the proc- 
ess of determining a curve in colorant space comprises the following steps : 

determining at least one 2-ink process from said printer model ; 

- determining at least one 1 -ink process from each determined 2-ink process ; 

in each determined 1 -ink process, determining at least one set of colorant values, within said colorant gamut ; 



15 



EP0763 930 A1 

. construing said curve in cdorant space, having two sets of colorant values. be.ong*9 to one 2-inK process, 
as endpoints. 

5 ised Neugebauer equations. 

ther comprising the steps of : 
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Fig. 6 
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